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Abstract 

Option pricing is mainly based on ideal market conditions which are well represented 
by the Geometric Brownian Motion (GBM) as market model. We study the effect 
of non-ideal market conditions on the price of the option. We focus our attention 
on two crucial aspects appearing in real markets: The influence of heavy tails and 
the effect of colored noise. We will see that both effects have opposite consequences 
on option pricing. 



1 Introduction 



Option pricing lies, from the very beginning, at the heart of mathematical fi- 
nance as scientific discipline. In 1900, Bachelier proposed the arithmetic Brow- 
nian motion for the dynamical evolution of stock prices as a first step towards 
obtaining a price for options [1]. Since then and specially after the landmark 
work of Black and Scholes [2] and Merton [3], the underlying asset of any 
option is usually assumed to be driven by Brownian motion, being the most 
ubiquitous model the geometric Brownian motion (GBM): 

§ = ^ + ^ 2 + (i) 



where S(t) is the price of stock (or the value of an index) at time t, fi+ o~ 2 /2 
is the drift, a is the volatility and £(t) is Gaussian white noise with zero mean 
and correlation (£(£)£(£')) = S(t — tf). Associated with stock price S(t) there is 
another random process: the return. In continuos-time finance this is denned 
by R(t) = In S(t)/S (to). For the GBM (1) one can see, after using Ito lemma 
[4], that the return obeys the following SDE 

R = fjL + a£{t), (2) 



Preprint submitted to Elsevier Preprint 



1 February 2008 



from which we see that [i is the mean return rate. Since, by definition, the 
initial return is zero, then in terms of the integral of £(£) which is taken in the 
Ito sense, we can explicitly write the return as 



where the integral of appearing in the rhs of this equation is simply the 
Wiener process W(t). Note that, independently of the stock price S(to), the 
return at initial time t is always zero since R(to) = In S (to) / S (to) = In 1 = 0. 

However, the market model given by Eq. (1) or Eqs. (2)-(3) does not fully 
account for all properties observed in real markets. Some of them are quite 
crucial such as the existence of heavy tails and negative skewness in probability 
distributions, evidence of self-scaling (at least for high frequency data) [5], 
leverage [6] and even the existence of mild correlations in the increments of 
return [7]. All of this, undoubtedly has to affect option pricing. Therefore, 
the theoretical price for the option obtained through the assumption that 
markets are governed by the GBM would not be adjusted to reality. Our main 
objective here is to elucidate how some of the above properties of real markets 
may affect the price of options. Specifically we will obtain an option price for 
market models that present fat tails and self-scaling. After that, we will relax 
the so called efficient market hypothesis [8] and assume the existence of auto 
correlations in stock prices. We restrict the analysis to European call options, 
although it can be easily generalized to the European put and other options 
with fixed exercise date. The extension to American options seems to be more 
involved. 

The starting point of this work is the martingale option price. As was shown 
in early eighties, the Black-Scholes (B-S) option price can also be obtained 
using martingale methods [9]. This is a shorter, although more abstract way, 
to derive an expression for the call price. The main advantage is that one only 
needs to know the probability density function governing market evolution. 
The main inconvenience is that there is no known hedging strategy as in the 
case of B-S theory. 

The probability distribution of some of the more realistic market models we 
will use is only known through its characteristic function. We will therefore 
derive a formula which will allow us to obtain, via martingale price, the op- 
tion price in terms of the characteristic function of the market model chosen. 
Having the necessary tools to address the problem, we will apply them to a 
recently presented market model which shows self-similarity, fat tails and fi- 
nite moments [10]. We will obtain an option price this model and compare it 
with B-S price. 




(3) 



to 
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The second problem we want to address is the effect of colored noise in the 
price of options. The GBM model supposes that prices are driven by Gaussian 
white noise, this implies that prices S(t), at different times, are uncorrected 
which, for certain times scales, does not seem to be completely adjusted to 
reality. We will thus assume that prices are driven by colored noise and obtain 
the corresponding fair price for the European call. 

The paper is organized as follows, after this introduction, in Sect. 2 we review, 
from a physicist's point of view, the martingale price for the European call. 
In Section 3 we obtain the call price formula using Fourier analysis and apply 
it in Sect. 4 to heavy tail and self-similar processes. In Sect. 5 we study the 
effect of colored noise on option pricing. Conclusions are drawn in Sect. 6. 



2 The martingale option price 

We start by recalling that an European call option is a derivative giving to its 
owner the right but not the obligation to buy a share at a fixed date T (the 
maturity) for a certain price K (the strike or striking price). The call option 
contract is basically specified by its gain at maturity. If S — S(T) is the share 
price at time T this gain, called payoff, is given by 

(S-K) + = max[S-K,0]. (4) 

As was shown in the early eighties [9] , the B-S option price can also be obtained 
using martingale methods. This is a shorter, although more abstract way, to 
derive an expression for the call price. The main advantage being that one only 
needs to know the probability density function governing market evolution. 

Before proceeding further let us briefly explain, in layman terms and without 
pretending to be rigorous, what a martingale is. Suppose X(t) is a well defined 
random process such that its previous history is perfectly known. That is, we 
keep record, and hence know of, all values taken by X(t) before time t. Under 
these conditions the process is a martingale if the conditional expected value 
{X(t)\X(t ) = x ) for any t > t is given by 

(X(t)\X(to)=x )=x , (t>f ). (5) 

An example. Under quite general conditions any driftless process driven by 
zero-mean noise is a martingale [11]. If X(t) follows the SDE: X = F(t) where 
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F(t) is a zero- mean random process, then 

t 

X(t) = x + J F{t')dt\ 

to 

where xo = X(t ) is known and F(t) is supposed to be integrable, at least in 
the sense of generalized functions. Since (F(t)) = 0, we have (X(t)\xo) = xq. 
Hence X(t) is a martingale. 

In what follows we generalize the GBM model (l)-(3) to include driving noises 
other than the Wiener process and thus trying to account for some features 
of real markets such as fat tails, self-scaling and, eventually, correlations. We 
therefore proposed as market model the one whose return is given by (see Eq. 
(3)) 

t 

R(t) = n{t - t ) + J F(t')dt', (6) 

to 

where \i is the (constant) mean return rate, F(t) is zero-mean noise and, if 
needed, the integral of F(t) should be interpreted in the sense of Ito. We define 
the zero-mean return X(t) by 

t 

X(t) = J F{t')dt'. (7) 

to 

In generalizing the GBM for option pricing we have to assume that the zero- 
mean return X(t) does not allow for arbitrage opportunities. In other words, it 
must avoid money profits for free without taking any risk. To this end, as was 
proved by Harrison et al. [12], it suffices that X(t) be a noise of unbounded 
variation. We will assume this condition in what follows. 

One can easily show from Eq. (7) that 

(X(t 2 )\X(t 1 )}=X(t 1 ) 

for any t 2 > t±. Hence, the zero-mean return is a martingale. Moreover, in 
terms of X(t) the return is 

R(t)=fi(t-t )+X(t), (8) 

in consequence, R(t) obeys the following SDE 

R(t)=fi + F(t), (9) 
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which is the generalization of SDE (2). 

We now return to options. The equivalent martingale measure theory imposes 
the condition that, in a "risk-neutral world", the stock price S(t) evolves, on 
average, as a riskless security and discards any underlying process that allows 
for arbitrage opportunities. Thus if at time t we buy a European call on stock 
having price S, then in terms of the risk-neutral density p*(R, t\to) it is possible 
to express the price for the European call option by defining its value as the 
discounted expected gain due to holding the call. That is [9] 

oo 

C(S,t) = e- riT - t) J (Se R - K)p*(R,T\t)dR, (10) 

where T and K are the maturity and strike of the option. 

Let us define the risk-neutral density p*(R,t\t ). This density, which is also 
called the "equivalent martingale measure" associated with return R(t) under 
"risk- neutrality" , is the probability density function (pdf) of a process R*(t) 
which we call risk-neutral return and such that 

(e**^) = e^-*°\ (11) 

where r is the risk-free interest rate ratio. Since the stock price associated with 
R* is given by S*(t) = S(to)e R * (:t \ condition (11) imposes that risk-neutral 
prices S* must grow on average as a riskless security. 

Note that we can easily prove that stock price S is always greater than call 
price. In effect, from (10) we have 

oo 

C(S,t) < Se- r V-V J e R p*(R, T\t)dR = Se~ r ^ (e R *^\t) 

— oo 

and, after using Eq. (11) we conclude that 

C{S,t)<S. (12) 

We now assume that the market model for the evolution of the underlying 
stock is given by Eq. (9). Then, analogously to Eq. (8) for the return R(t), we 
impose that the risk-neutral return R*(t) be written as 

R*(t)=m*(t-t )+X(t), (13) 
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where m*(t) is a deterministic function and X(t) is given by Eq. (7). In order 
to find the expression for m*(t), observe that 



e R*(t)\ = e m*(t) / e x(t) 



and from the risk neutrality condition (11) we see that 

•(t-t )=r(t-t )-\n(e x ^). 



m 



The average (e x ^) is obtained from the known probability distribution of 
the driving noise X(t). Indeed, let (px(u,t\t ) be the characteristic function 
of X(t), 

Vx (uj,t\t )^(e^\t ), (14) 



then 



e x(t) ) = <px(-i,t\t ) (15) 



and 



m*(t-t )=r(t-t )-kiip x (-i,t\t ). (16) 



For the GBM, the zero- mean return X(t) is the Wiener process (cf. Eq. (2)). 
Thus 

Vx (u,t\t ) = e'^ t - t ^ 2 , (17) 



consequently In (px(— h t\to) — °" 2 (^ — to)/2 is the spurious drift, and 

m*(t-t ) = (r-a 2 /2)(t-t ). (18) 

Finally, since R*(t) and X(i) are related by the simple linear relation (13), 
the risk-neutral density p* can be written, in terms of the zero-mean return 
density p x , as 

p*{R,t\t ) =Px{R-m*{t-t ),t\t ), (19) 
and the price of the call is given by 

oo 

C(S, t) = e - r ( T -*> J (Se R - K)p x (R-m*(T -t),T\t)dR. (20) 

\n(K/S) 
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3 Martingale option price by Fourier analysis 

As we have said, the main advantage of martingale pricing is that one can 
obtain a fair option price when the market obeys a random dynamics different 
than the geometric Brownian motion. In this case one only has to know the 
return distribution of the underlying and then, using Eq. (20), one readily gets 
a fair price for the option. Nevertheless, knowing an analytical expression of 
the pdf px(x,t\t ) may be, in practice, beyond our reach. There are however 
many situations in which we know the characteristic function of the model 
although one is not able to invert the Fourier transform and obtain the density 
Px(%, t\to). This is the case of non-Gaussian models such as Levy processes [5] 
and some generalizations [10], the case of stochastic volatility models [13,14] 
among others. For these cases, we will develop an option pricing based on a 
combination of martingale methods and harmonic analysis. 

We recall that the characteristic function (cf) of a random variable is the 
Fourier transform of its probability density function. Thus the characteristic 
function of the return will be given by: 



The starting point of our analysis is Eq. (10) that we write in the form 



The first integral on the right hand side (rhs) of this equation can be solved 
in closed form with the result (see Eq. (11)) 



oo 




— oo 




(21) 



oo 




— oo 



As to the second integral on the rhs of Eq. (21), we have 



\n(K/S) 



oo 




(Se R - K)p*{R, T\t)dR = K 



j (e- z - l)p*(-z - R K ,T\t)dR, 



— oo 
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where 



R K = ln(S/K) 



(22) 



can be considered as the return associated with moneynes^] The inverse 
Fourier transform of the characteristic function 



p *(R,T\t) = ^- I e^ R ^{u : T\t)du 



allows us to write this second integral as 



I=Y I V*(u,T\t)e iulRK duj J e iu)Z (e~ z - l)dz. 

-oo 



(23) 



In the Appendix A we show that / can be written as 



17 + 9~~ / <P (u,T\t)- ; I 

2 2ir J 1 — iuj 

— oo 



e wHK <p*(u,T\t) - 1 —-(24) 

J L) 



duo 



Collecting results into Eq. (21), we finally get 



C(S,t) = S- |e-^) - f J <p*( Ut T\t)^ 



RK d 



to 



—i 



e iujRK (p*(uj,T\t) - 1 



- IUJ 

duj 1 
oo I 



(25) 



The representation (25) is indeed very useful when the risk-neutral density p* 
is unknown but its characteristic function (p* is known. This would be indeed 
the case of more sophisticated market models such as those of next sections. 
A similar result has been given in [15] using a different but equivalent form of 
Eq. (25) in order to numerically perform the integration knowing (p*(u,T\t). 
It is asserted in [15] that these Fourier methods allow a fast computing of the 
option price. 

We close this section presenting an alternative price formula to that of Eq. 
(25) that is written in terms of the characteristic function of the zero-mean 
return instead of the risk- neutral cf ip*. This price formula only involves real 



1 The neologism "moneyness" refers to the ratio S/K. 
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quantities and it can be more convenient in a number of cases. From Eq. (19) 
we easily see that 

ip*{u,T\t) = e-'^V^.Tjt), (26) 

where ipx(oo,T\t) is the characteristic function of the zero-mean return X(t). 
Moreover, since (X(t)\t ) = 0, the distribution of X(t), px(x,T\t), is symmet- 
ric around x = 0. Then its Fourier transform is a real and even function of uj 
[17], 

<p x {-u,T\t) = ip x {u,T\t). 

These properties allow us to further simplify the integrals appearing in Eq. (25). 
In effect, in terms of tpx(w, T\t) given by Eq. (26), the first integral on the rhs 
of Eq. (25) (we call it I±) reads 

n °f , rTnl . cos coa(T — t) — to sin cua(T — t) , 
h = 2 / <p x (w,T\t) ^ > —— 2 ^ '-duj, 

J 1 + LO Z 



where 

a(t) = R K + m*(t). (27) 



Following an analogous reasoning we obtain, for the second integral on the rhs 
of Eq. (25) the result 

oo . . 

f / n-i i \ sin ua(T — t) , 
I 2 = 2i <px(w,T\t) ^ '-duo. 

J LU 


The substitution of these two integrals into Eq. (25) yields our final result 



C(S,t) = S-^e 



-r(T-t) 



K 



7T 



-r(T-t) 



<p x (w,T\t) 



cos ua(T — t) + 



sina;a(T — t) 



du 



(28) 



where a is given by Eq. (27). Note that this expression only involves one real 
integral and it is therefore simpler and faster to compute than Eq. (25). 

Finally, for the GBM market model, the substitution of Eqs. (17)— (18) and 
(27) into Eq. (28) yields the Black-Scholes formula 



C BS (S,t) = SNidJ - Ke- rt N(d 2 ), 



(29) 
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where 

X 

N{x) = (l/v 7 ^) J e~ z2/2 dz } 

— oo 

is the probability integral, and 

rt + R K + a 2 /2 



do 



rt + R K - a 2 /2 



(30) 
(31) 



In Fig. 1, we plot this classic option price using the volatility estimated from 
Standard and Poors-500 tick by tick data (January 1986 - December 1996). 
Time to maturity is assumed to be 30 days and we compare the resulting call 
price with the deterministic call price given by Eq. (55). 
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Fig. 1. The B-S call price in terms of the moneyness. For this graph, we take 
r = 5% year -1 and the volatility estimated from the S&P-500. The B-S call price 
formula is given by Eq. (29) and the deterministic price is given by Eq. (55). 



4 Option pricing on heavy-tailed stocks 



We will now apply the formalism developed hitherto to pricing options whose 
underlying assets have probability distributions showing fat tails. In such 
cases, the ideal market conditions upon which the B-S theory is based are 
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no longer valid and B-S formula is not reliable. We want to elucidate what 
this implies on the price of the option. 



The first market model incorporating heavy tails in a natural way was given 
by Mandelbroth in the early sixties [18,19] who, based on Pareto-Levy stable 
laws, obtained a Leptokurtic distribution. In that model, the zero-mean return 
defined above has the following characteristic function 



(1 < a < 2, k > 0). There is, however, a drawback in this model: no finite mo- 
ments exist beyond the first one. This certainly is a severe limitation specially 
for the practitioner. Moreover, testing the Pareto-Levy distribution against 
data has always resulted with the same conclusion: the tails are far too heavy 
compared with actual data. Although, as Mantegna and Stanley have clearly 
shown [5] , the Levy distribution fits very well the center of empirical distribu- 
tions - much better than the Gaussian distribution of the GBM. Furthermore, 
Levy distributions show a feature recently observed in data, specially in intra- 
day data, and this is the self-scaling behavior of the probability distribution at 
different times [5,20,22]. Nevertheless, and regardless these advantages, the ab- 
sence of finite moments seriously affects obtaining a price for options. Indeed, 
in this case the average {e x ^) does not exist (see Eqs. (15)— (16)) and option 
price (20) is useless. Therefore, any attempt to price options through Levy 
processes has to be done by means of truncated Levy distributions [24] with 
the added inconvenience that truncation is an ad hoc procedure introducing 
high degrees of arbitrariness. 

In a recent work we have presented a market model that explains the appear- 
ance of fat tails and self-scaling but still keeping all moments finite [10,23]. In 
that model we assume that all random changes in the stock price are modelled 
by a continuos superposition of different shot-noise sources, each source cor- 
responding to the detailed arrival of information [10]. The zero-mean return 
is given by 



where n(t) represents the maximum "number" of noise sources taken up to 
time t and it is an increasing function of time to be determined by the self- 
scaling property (see below). For any fixed time t, Y(u,t) are independent 
random variables for different values of parameter u and for any fixed value of 
t, Y(u, t) is a white shot-noise process represented by a countable superposition 



ipx(uj,t) = exp(—ktu a ), 



(32) 



n(t) 




(33) 



— oo 
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of Poisson pulses of rectangular shape: 

oo 

Y(u,t) = Y i A k ( U )e[t-T k ( U )], (34) 



k=i 



where 6(t) is the Heaviside step function, T k (u) marks the onset of the fcth 
pulse, and A k {u) is its amplitude. Both T k {u) and A k {u) are independent and 
identically distributed random variables with probability density functions 
given by h(x,u) and ip(u,t), respectively. 

We assume that the occurrence of jumps is a Poisson process, in this case 
the shot-noise Y(u, t) is Markovian, and the pdf for the time interval between 
jumps, ifj(u,t)dt = P{t < T k {u) — T k _i{u) < t + dt}, is exponential: 

i/>(u,t) = A(w)e- A ^' (t>0), 

where X(u) is the mean jump frequency. We recall that jump amplitudes A k (u) 
are identically distributed (for all k = 1, 2, 3, • • •) and independent random 
variables (for all k and u). In what follows we will assume that they have zero 
mean and a pdf, h{u,x)dx = P{x < A k (u) < x + dx}, depending on a single 
"dimensional" parameter which, without loss of generality, we assume to be 
the standard deviation of jumps a(u) = y {A\{u)). That is, 



h(u, x) = h 
a(u) 



x 



a(u) 



The function n(t) appearing in Eq. (33) is determined by imposing on X(t) 
self-scaling properties. It reads (see [10] for details) 

n(t) = -\nbt. (35) 
a 



Under these assumptions we have shown that the characteristic function of 
the zero mean return X(t) is [10] 



tpx(u,t) = exp < 



-abt 



z- L - a [l - h(ujz)]dz 



(36) 



where h(cu) is the characteristic function of the jump distribution h(x). We 
thus see that the model depends on three positive constants: a, b, and a < 2 
and an unknown distribution h(x) to be guessed from data. 

Let us summarize the main results and consequences of the model. First, the 
model contains the Levy process as a special case. Indeed if a — > and b — > oo 
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in such a way that ah is finite, Eq. (36) yields, for any jump distribution h(x), 
the Pareto-Levy distribution (32) with 

oo 

k = ab J z~ l - a [l - h{z))dz. 
o 

Since the limit b — > oo implies that n(t) — > oo (cf. Eq. (35)) we see from Eq. 
(33) that Levy process are continuos an unbounded superposition of families 
of white Poissonian shot noises^. 

Second, the volatility of the zero-mean return is given by 

(X 2 (i)> = -i- {btfl^ (37) 
2 — a 



which proves that a < 2 and that the volatility shows super-diffusion. The 
anomalous diffusion behavior of the empirical data (at least at small time 
scales) was first shown by Mantegna and Stanley [25] . Third, kurtosis is con- 
stant and given by 

^ = (2 a)^»(0) 
(4 — a)a 



Thus 72 > for all t, in other words, we have a Leptokurtic distribution (i.e., 
heavy tails) in all time scales. Fourth, the return probability distribution scales 
as 

<Px(u,t)=<px((bt)V a u) (39) 



and the model becomes self-similar [5,20,22]. 

As to the asymptotic behavior of our distribution. It can be shown from 
Eq. (36) that the center of the distribution, defined by \x\ < (bt) 1 ^, is again 
approximated by the Levy distribution defined above. On the other hand the 
tails of the distribution are solely determined by the jump pdf h(u) by means 
of the expression 

oo 

/ U "H^ (M » ( 6t ) 1A *)- ( 4 °) 

\x\/a m 



2 We use the adjective "unbounded" in the sense that there is no maximum number 
of shot-noise sources, i.e., n(t) = oo in Eq. (33). 
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-15 -10 -5 5 10 15 

x/L 

Fig. 2. Probability density function p(x, t) for t = 1 min. Circles represent empirical 
data from S&P 500 cash index (January 1988 to December 1996). £ is the standard 
deviation of empirical data. Dotted line corresponds to the Guassian density. Dashed 
line is the Levy distribution and the solid line is the Fourier inversion of Eq. (36) 
with a gamma distribution of jumps. 

Therefore, return distributions present fat tails and have finite moments if 
jump distributions behave in the same way. This, in turn, allows us to make 
statistical hypothesis on the form of h{u) based on the empirical form and 
moments of the pdf. 

In Fig. 2, we plot the probability density px{x,t) of the S&P 500 cash index 
returns X(t) observed at time t — 1 min (circles). S = 1.87 x 10~ 4 is the 
standard deviation of the empirical data. Dotted line corresponds to a Gaus- 
sian density with standard deviation given by E. Solid line shows the Fourier 
inversion of Eq. (36) with a = 1.30, a m = 9.07 x 10" 4 , and a = 2.97 x 10~ 3 . 
We use the gamma distribution of the absolute value of jump amplitudes, 

h(u) = ^\uf- 1 e-^/2T{(3), (41) 



with p = 2.39, and /x = J (3(/3 + 1) = 2.85. Dashed line represents a sym- 
metrical Levy stable distribution of index a = 1.30 and the scale factor 
k = 4.31 x 10~ 6 [10]. We note that the values of a m and S predict that 
the Pareto-Levy distribution fails to be an accurate description of the empir- 
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ical pdf for x 3> 5S (see Eq. (40)). We have chosen a gamma distribution of 
jumps because (i) as suggested by the empirical data analyzed, the tails of 
p(x,t) decay exponentially, and (ii) one does not favor too small size jumps, 
i.e., those jumps with almost zero amplitudes. 

We will now obtain a price for a European call assuming that the underlying 
has a zero-mean return given by Eq. (33). In order to apply the formalism 
derived in Sect. 2 and 3 we have to see whether our X(t) allows arbitrage, in 
other words, we have to prove that X(t) has unbounded variations (see Sec. 
2 and Ref. [12]). We proceed as follows, the variation of X(t) is given by its 
time derivative X(t). From Eqs. (33)-(34) we have 



n(t) 



X(t)=n{t)Y(n(t),t)+ / A k (u)5[t-T k (u)}. 



Following an analogous procedure to that of Refs. [10] and [23] we can see 
that the auto correlation of the derivative process is, in terms of the volatility 
(X 2 (t)> (cf. Eq. (37)), given by 



(X(t + r)X(t)) 



(X 2 (t)} 



(2 — a) a 



at 



+ S(t) 



which clearly shows that X(t) is a process of unbounded variation. 

Assuming a gamma distribution of jumps, Eq. (41), it can be proved that the 
characteristic function of X(t) given by Eq. (36) explicitly reads 



(px(w,t) = exp 



a 



1 — Re < F I /3, —a; 1 — a; 



iuj{bt) 



l/a N 



(42) 



where F(a,b; c; z) is the hypergeometric function and Re(-) is the real part. 
The substitution of this equation into Eq. (28) results in the price of the call. 
Unfortunately, one cannot evaluate analytically the resulting integrals and 
numerical work should be performed. 

In Fig. 3, we plot the call price CT a u{S,t) obtained by numerical integration 
of Eq. (28) with Eq. (42) in terms of moneyness S/ K. We compare this price 
with B-S price with the same time to maturity, T — t, volatility and risk- 
free interest rate. We see that the inclusion of heavy tails in the probability 
distribution substantially increases option prices. In Fig. 4, we plot CT a u(S,t) 
as a function of S/ K for three different times to maturity. We confirm that in 
this case (as in B-S case) an increase of time to maturity implies an increase 
of the call price. 
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Fig. 3. The heavy tail and B-S call prices in terms of the moneyness when maturity 
is in 30 days. For this graph, we take r = 5% year -1 and the parameters estimated 
for the S&P-500 from January 1988 to December 1996 (a = 3.54 x lCT 3 days~ 1/2 ). 
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Fig. 4. The heavy tail call price in terms of the moneyness. For this graph, we take 
r = 5% year -1 and the parameters estimated for the S&P-500. 
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Fig. 5. The relative differences, A = {Cxaii — C bs) I CTaiii m terms of the moneyness 
S/K (r = 5% year' 1 , a = 3.54 x 10- 3 days~ 1/2 ). 

Finally in Fig. 5 we plot the relative difference 

A = {Cxail ~ CBs)/CTail 

as a function of moneyness. We observe that for OTM, ATM, and some ITM 
options A increases as T — t increases, while for options well in the money A 
becomes smaller with T—t. There is thus a "crossover" which, in our example, 
is located around S/K = 1.03. We also note that A — > as S/K —>■ oo very 
slowly which implies the great persistence of the heavy-tail effect in the price 
of the option. 



5 Option pricing on correlated stocks 

As we have mentioned in Section 1, the classical option formula of Black, Sc- 
holes and Merton is based on the "efficient market hypothesis" which means 
that the market incorporates instantaneously any information concerning fu- 
ture market evolution [8]. However, as empirical evidence shows, real markets 
are not efficient, at least at short times [26]. Indeed, market efficiency is closely 
related to the assumption of totally uncorrelated price variations (white noise). 
But white noise is only an idealization since, in practice, no actual random 
process is completely white. For this reason, white processes are convenient 
mathematical objects valid only when the observation time is much larger than 



17 



the auto correlation time of the process. And, analogously, the efficient mar- 
ket hypothesis is again a convenient assumption when the observation time is 
much larger than time spans in which "inefficiencies" (i.e., correlations, delays, 
etc.) occur. 

Note that auto correlation in the underlying driving noise is closely related to 
the predictability of asset returns, of which there seems to be ample evidence 
[28,29]. Indeed, if for some particular stock the price variations are correlated 
during some time r, then the price at time t 2 will be related to the price at a 
previous time t\ as long as the time span t 2 — t\ is not too long compared to 
the correlation time r. Hence correlation implies partial predictability. 

In Ref. [27], we have derived a nontrivial option price by relaxing the effi- 
cient market hypothesis and allowing for a finite, non-zero, correlation time 
of the underlying noise process. As a model for the evolution of the market 
we chose the Ornstein-Uhlenbeck (O-U) process for essentially two: (a) O-U 
noise is still a Gaussian random process with an arbitrary correlation time 
r and it has the property that when r = the process becomes Gaussian 
white noise, as in the original Black-Scholes option case, (b) The O-U process 
is, by virtue of Doob's theorem, the only Gaussian random process which is 
simultaneously Markovian and stationary. In this sense the O-U process is the 
simplest generalization of Gaussian white noise. 

We thus assume that the underlying asset price is not driven by Gaussian 
white noise but by O-U noise V(t). In other words, instead of Eq. (2), we 
now have the following singular two-dimensional diffusion process: 

R(t)=» + V(t), (43) 
V(t) = \ l-V(t) + aft)] , (44) 

where r > is the correlation time. More precisely, V(t) is O-U noise in 
the stationary regime, which is a Gaussian colored noise with zero mean and 
correlation function: 

(VihMh)) = ^- e -\^\/\ (45) 
It 

Note that, when r = 0, this correlation goes to <J 2 5(t 1 —t 2 ) and we thus recover 
the one-dimensional diffusion discussed above. Therefore, the case of positive 
r is a measure of the inefficiencies of the market. The combination of Eqs. (43) 
and (44) leads to a second-order stochastic differential equation for R(t) 

TR(t)+R(t)=n + <r£(t). (46) 
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From this equation, we clearly see that when r = we recover the one- 
dimensional diffusion case (2). In the opposite case when r = oo, Eq. (44) 
shows that V(t) = 0. Thus, V(t) is constant which we may equal it to zero. 
Hence, R(t) = fit and the underlying price, S(t) = Soe^, evolves as a riskless 
security. We also observe that the O-U process V(t) is the random part of the 
return velocity, and sometimes we refer to V(t) as the "velocity" of the return 
process R(t). 

One may argue that the O-U process (43)-(44) is an inadequate asset model 
since the return R(t) given by Eq. (43) is a continuous random process with 
bounded variations. As we have explained in Section 2, continuous processes 
with bounded variations allow arbitrage opportunities and this is an undesir- 
able feature for obtaining a fair price. Thus, for instance, in this present case 
arbitrage would be possible within a portfolio containing bonds and stock 
whose strategy at time t is buying (or selling) stock shares when fi + V{t) is 
greater (or lower) than the risk-free bond rate [12]. However, the problem is 
that in the practice the return velocity V(t) is non tradable and its evolution 
is ignored. In other words, in real markets the observed asset dynamics does 
not show any trace of the velocity variable. Indeed, knowing V(t) would imply 
knowing the value of the return R(t) at two different times, since 

rr/ , R(t)-R(t-e) 
V(t)= hm -±> ^ — fx. 

e^0+ e 

Obviously, this operation is not performed by traders who only manage port- 
folios at time t based on prices at t and not at any earlier time. This feature 
allows us to perform a projection of the two-dimensional diffusion process 
[R(t),V(t)] onto a one-dimensional equivalent process R(t) independent of 
the velocity V(t). We have shown elsewhere [27] the projected process R(t) is 
equal to the actual return R(t) in mean square sense, ([R(t) — R(t)} 2 ) = for 
any time t, and it obeys the following one-dimensional SDE[]] 

R(t) = »+y/k(T-t)Z(t), (47) 

where 

k{t) = a 2 (l- e-*/ T ) . (48) 

In this way, we have projected the two-dimensional O-U process (S, V) onto 
a one-dimensional price process which is a Wiener process with time varying 
volatility. Therefore, the return given by Eq. (47) is driven by a noise of un- 
bounded variation, the Wiener process, and the O-U projected process is still 

3 Since R(t) and R(t) are equal in mean square sense we will drop the bar on R as 
long as there is no confusion. Thus, we will use R for the projected process as well. 
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a suitable starting point for option pricing since it does not permit arbitrage. 
We also note that we need to specify the final condition of the process because 
the volatility \fk is a function of the time to maturity T — t, and this implies 
that the projected asset model depends on each particular contract. 

The zero-mean return (7) associated with the projected process (47) follows 
the SDE 



X(t) = y/k(T-t)t(t), 
and its density, px(x,T\t), satisfies the backward equation 

dpx l. frp ,x d 2 p x 

with final condition p x (x,T\T) = S(x). The solution to this problem reads 

1 



p x (x,T\t) 



^2hk(T - 1) 



cxp 



X 



2k(T - 1) 



(49) 



where (cf. Eq. (48)) 



K (t) = a 2 [t - r (l - e*/ r ) 



(50) 



The corresponding characteristic function reads 

<Px(w,T\t) = exp[-/c(T - t)u 2 /2], 
and from Eq. (16) we get 

m*(T -t)= r(T -t)- k(T - t)/2. 



(51) 



In this case the price of the European option is obtained substituting Eqs. 
(49) and (51) into Eq. (20). We have 



C ou (S,t) = S N(d° u ) - Ke- r ^ N(d° u ), 



(52) 



where N(z) is the probability integral, and 



d o V = HS/K) + r(T - t) + k(T - t)/2 



Ik(T - 1) 



(53) 
(54) 



with n(t) given by Eq. (50). 



20 




moneyness 

Fig. 6. The correlated and B-S call prices in terms of the moneyness when maturity 
is in 15 days and correlation time is 2 days. For this graph, we take r = 5% year -1 
and the parameters estimated for the S&P-500 (er = 3.54 x lO^days -1 / 2 ). 

When r = 0, the variance becomes n(t) = a 2 t and the price in Eq. (52) reduces 
to the Black-Scholes price given by Eq. (29). We observe that the O-U price 
in Eq. (52) has the same functional form as B-S price in Eq. (29) when a 2 t is 
replaced by n(t). 

In the opposite case, r = oo, where there is no random noise but a determin- 
istic and constant driving force (in our case it is zero), Eq. (52) reduces to the 
deterministic price (cf. Eqs. (4) and (10)) 

C d (S, t) = max[S - Ke~ riT - t} , 0]. (55) 



We have proved elsewhere [27] that the O-U price is an intermediate price 
between B-S price and the deterministic price (see Fig. 6) 

C d (S,t)<C ou (S,t)<C BS (S,t), (56) 



for all S and < t < T. 

In Fig. 7, we plot Cou m terms of S/ K for three different values of the time to 
maturity. We see in the figure that, analogously to B-S and heavy tail cases, 
the O-U price increases with maturity. In any case, the assumption of uncor- 
related underlying assets (B-S case) overprices any call option. This confirms 



21 



- 1 1 

correlated call -maturity 1 5 d.- 
correlated call -maturity 30 d.- 
correlated call -maturity 60 d.- 




0.96 0.98 1 1.02 1.04 1.06 1.08 

moneyness 

Fig. 7. The correlated call prices in terms of the moneyness with correlation time 
of 2 days. For this graph, we take r = 5% year^ 1 and the parameters estimated for 
the S&P-500 (a = 3.54 x l(T 3 days" 1/2 ). 




moneyness 

Fig. 8. Relative differences, D = (Cbs — Cou)/Cbs-, m terms of moneyness (r = 2 
days, r = 5% year -1 , and a = 3.54 x 10~ 3 days -1 / 2 ). 
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the intuition understanding that correlation implies more predictability and 
therefore less risk and, finally, a lower price for the option. 

In fact, one can quantify this overprice by evaluating the relative difference 

D = (Cbs — Cou)/Cbs- 

Figure 8 shows the ratio D(S,t), for a fixed time to expiration, plotted as a 
function of the moneyness, S/ K. Note that the relative difference D increases 
as T—t decreases. Moreover, contrary to the fat tail case, D quickly fades away 
with moneyness and it is negligible for ITM options showing no crossover. 



6 Conclusions 

We have performed an analytical study of option pricing on a more realistic 
setting than that of Black and Scholes theory. Our starting point has been 
the martingale option price conveniently modified to become more practical to 
implement by means of Fourier analysis. We have applied these analytical tools 
to study the effect on option price of long tails in the underlying probability 
distribution. For this, we have used a market model we recently developed 
[10] which combines fat tails, finite moments and self-similarity. The second 
problem we wanted to address is the effect of colored noise on option pricing, 
since there seems to be empirical evidence of mild auto correlation in stocks 
[7]. We have thus rederived, using the methods herein presented, our previous 
result on option pricing on stocks driven by O-U noise. 

One overall conclusion is that, taking the Black-Scholes price as a benchmark, 
heavy tails overprice by far options, while colored driving noise underprices it. 
The scheme is 

S > C Ta ii > Cbs > C u > C d , 

where S is the underlying price when the option is bought and Cd is the 
deterministic price given by Eq. (55). Moreover, the relative difference with 
B-S price is in both cases more important for OTM and ATM options than 
for ITM options. In any case, there is a great persistence of fat-tail effects on 
the price of the options. 

There are some open questions, one of which being the extension of the above 
formalism to other option classes (American, Asian, etc.). Another point is to 
investigate whether fat tails and correlations produce the "smile effect" in the 
implied volatility. We believe that this is more likely to happen in the case of 
heavy-tails. This point is under present investigation. 
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A Derivation of Eq. (24) 



Let us derive Eq. (24). We start from Eq. (23) 



K 
2^ 



oo oo 

J ip*{uj,T\t)e iujRK duj J e iu}Z (e~ z -l)dz, 



but 



and 



OO CXJ 

J e iwz (e- z -l)dz = — — - J e iuJz dz, 
o o 



oo oo 

/ e iuJz dz = lim / e'^'^dz = lim 

J e^0+ J e->(H 



o+ to + is lo + iO ' 
u u 

where this result has to be understood in the sense of generalized functions as 
[16] 



= =R7r5(cj) + P 



where 5{oj) is the Dirac delta function, and P[l/u;] is the Cauchy principal 
value, i.e., 

P[l/u] = 1/u for u ^ 

and 



P 



duo 



lim 



-doj + 



-du 



LO 



where (f>(u) is any regular function fast decaying at infinity. Therefore, 

oo 

J e wz dR = t:5(u) + iP 



Hence 



e iwRK du 



RK <P*(u;,T\t) 



du 



u 



24 



and using the following alternative expression for the Cauchy principal value 
[16] 



cf>{u) 



du 



U 



</>H - 0(0) 



du, 



we obtain Eq. (24): 



RK du 



du 
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